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ABSTRACT 

We present the discovery of the widest (^-^ 6700 AU) very low mass field binary to date, found in a proper 
motion cross-match of the Sloan Digital Sky Survey and the Two Micron All Sky Survey. Our follow-up J-band 
imaging provides a 10-year baseline for measuring proper motions. Consequently, we are able to confirm the 
common proper motion of the pair to within 10 mas, implying a 99.5% probability of their physical association. 
Near infrared spectra of the components indicate spectral types of M6±l and M7±l. The system resides 
at a spectroscopic distance of 105 ± 13 pc and has an angular separation of 63.38 ±0.05". We have used 
evolutionary models to infer component masses of 0.105;!;[5y[7 Mq and 0.091;!;jjQ07 Mq. The large separation 
and low binding energy of this system can provide constraints for formation models of very low mass stars. 

Subject headings: binaries: general — stars: formation — stars: individual (2MASS J12583501H-4013083, 
2MASS J12583798H-4014017) — stars: low-mass, brown dwarfs 
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1. INTRODUCTION 

Studies of the low end of the stellar initial mass function 
have made great strides in recent years, and have revealed a 
large population of very low mass (VLM, M <0.1 Mq) stars 
and brown dwarfs (BDs) in the Galaxy. Their ubiquity poses 
a challenge to traditional models of star formation via grav- 
itational fragmentation, because their masses are well below 
the typical Jeans mass in molecular clouds. Whether these 
cool objects form similarly to more massive stars or require 
additional processes is intensely debated (e.g.lBureasser et alJ 
I2007t iLuhman et all 120071: IWhitworth et all l2007h . Since 
multiplicity properties can provide insight into formation sce- 
narios, it is of considerable interest to determine whether the 
properties of VLM binaries differ from those of more mas- 
sive stars. T rends pertaining to VL M binaries are discussed 
at length by iBurgasser et al.l (l2007h . These systems are ob- 
served to be less frequent, more tightly bound, and of higher 
mass ratios than their more massive counterparts. It is not 
clear whether these properties form part of a continuous mass 
dependent trend, or represent a unique population with a po- 
tentially different formation mechanism. 

In recent years much attention has been paid to relatively 
rare wide (> 100 AU) VLM binaries since their large sep- 
arations and low binding energies provide direct constraints 
for formation models. In particular, their existence challenges 
the ejection model for the formation of VLM stars and BDs 
(iReipurth & ClarkdllOOlh since such fragile systems are not 
expected to survive the ejection process. 

To date there are 8 known widely separated VLM and 
BD binaries, b o th in low-density star forming re gions (e.g. 
LuhrnanI 120041: iJava wardhana & Ivanov 2006; AUers et alJ 
20061 IClose et al.' '2007; Luhman et al i2009i) and in the 
field (e.g. Billeres et al. 2005; Caballer oH2007aHArtigau et alj 
I2008t FOizis et al.1 GOOOc iBejar et al.1 12008h . Surveys indi- 
cate a wide compan ion fraction for VLM stars and BDs of 
no more than 1-2% (IBurgasser et alj|2007l; ICaballeroll2007bl; 
I Allen & Reidll2008l) . However, searches for ultracool bina- 
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ries at the largest separations remain incomplet e. The re- 
cent d iscovery of the widest VLM binary by Artigau et alj 
( l2007h with a separati on of 5 100 AU, and the discovery of a 
~ 1 800 AU binary by lCaballerol (l2007ah point to the existence 
of a population of such ultrawide systems. 

Here we report the discovery of the widest VLM binary 
system to date with a projected separation of ~ 6700 AU. 

2. DISCOVERY AND OBSERVATIONS 

The system, comprised of 2MASS J12583501H-4013083 
and 2MASS J12583798H-4014017 (2M1258AB hereafter), 
was found in a cross-match of the Sloan Digital Sky 
Survey (SDSS) 6th data rele ase Photoprimary Catalog 
(lAdelm an-McCarthv et al.'2008) and th e Two Micron All Sk 
Survey (2MASS) Point Source Catalog dSkrutskie et al. l200 ^ 
in which we searched for common proper motion pairs con- 
taining VLM components. The correlation of catalogs, calcu- 
lation of proper motions, and the identification of co-moving 
stars was done in overlapping sections of 4 deg^ of the sky 
at a time, spanning the entire contiguous region of the SDSS 
Legacy survey in the northern galactic cap. For every 2MASS 
source the closest SDSS match was found and proper motion 
vectors with uncertainties were computed. A cut was made in 
order to select only stars that had moved at the 3(7 level com- 
pared to all other stars within the area. Stars within 120" of 
one another with proper motion amplitudes agreeing within 
2(7 and proper motion components agreeing within 1 cr in one 
of right ascension and declination, were flagged as potential 
binaries. The pair 2M1258AB was found by applying fur- 
ther color cuts of z' - 7 > 1.5, corresponding approximately to 
mid-M spectral types and later, to both components of poten- 
tial binaries within our subsample. 

To determine spectral types for the components of 
2M1258AB we obtained near-infrared spectra (R^250) on 
2008 March 01 using the SpeX Medium-Resolution Spectro- 
graph (Rav ner et alJ i2003) at NASA's Infrared Telescope Fa- 
cility (IRTF). Observations were made in prism mode with the 
0.3" sUt. We obtained six 120 s exposures arranged in three 
AB cycles. For telluric and instrumental transmission correc- 
tion, the AO star HD 109615 was observed immediately after 
the target at the same airmass. The spectra \yere re duced using 
SpeXtool ( Cushing et al]l200irv^cca"et alJl2003h . 

To confirm the common proper motion of this system. 
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TABLE 1 

Multiple EPOCHS of imaging observations for 2M1258AB 



Source 



Epoch (yr) 



a A (deg) 



Q B (deg) 



5 A (deg) 



<5 B (deg) 



2MASS 
SDSS 
CPAPIR 
CPAPIR 



1998.2666 
2003.3149 
2008.2292 
2008.4999 



194.645879(0.06) 
194.645997(0.05) 
194.646157(0.08) 
194.646156(0.08) 



40.218994(0.08) 
40.218824(0.04) 
40.218670(0.08) 
40.218668(0.08) 



194.658278(0.06) 
194.658460(0.06) 
194.658563(0.08) 
194.658549(0.08) 



40.233826(0.08) 
40.233639(0.06) 
40.233519(0.08) 
40.233522(0.08) 



Note. — SDSS and 2MASS positions are taken directly from tlieir respective catalogs. The uncertainties are given in tiie brackets, in 
units of arcseconds. 



follow-up y-band imaging observations were obtained on 
2008 March 25 and 2008 July 01 with the wide field near- 
infrared camera CPAPIR (E. Artigau et al., in preparation) at 
the 1.6-m Mont-Megantic telescope. At each visit, 16 expo- 
sures of 8.12 s were obtained by dithering the telescope by 
^5" between exposures. The follow-up /-band imaging data 
were reduced as follows. A sky frame was first constructed 
by taking the median of all images after masking the stars in 
each image. After subtraction of this sky frame, the images 
were divided by a normalized flat field image. The distor- 
tion and astrometry solutions of each image were computed 
by using the 2MASS PSC as a reference frame. The reduced 
and calibrated images were then combined by taking their me- 
dian. Astrometric uncertainties of ^ 80 mas were determined 
from the dispersion of the positions of all stars in the field 
with respect to their positions in the 2MASS PSC. The SDSS 
astrometry used was taken directly from the SDSS catalog. 
The positions reported by the catalog are derived from the r' 
band in which our sources are very faint. Therefore, as an 
additional check we calculated SDSS positions directly from 
the /' and z' images, and found them to be consistent with the 
catalog positions within the quoted uncertainties. A list of 
all astrometric measurements from OMM, SDSS and 2MASS 
used to compute proper motions are listed in table [T] 

Proper motions and their uncertainties were obtained using 
an error-weighted linear regression, yielding proper motions 
for the primary and secondary respectively of 76 ± 8 mas yr"' 
and 76 ± 8 mas yr"' in right ascension, and -1 15 ± 9 mas yr"' 
and -108 ± 10 mas yr"' in declination. The proper motions 
are in agreement within the Icr errors. The mean angular sep- 
aration of the system was found to be 63.38 ± 0.05". 

Using the s pace density for M dwarfs given by 
iPhan-Bao et alj ( l2003h and our total search area of 7668 
deg^, we would expect to find 1.8 pairs of M6-M8 dwarfs 
within 64" of one another and at distances ranging from 75- 
125 pc (spectroscopic distances are derived in |3j within 
our search area. Using the space velocity distributions from 
iBochanski et al.l (l2005l) the probability of any two M dwarfs 
having proper motions greater than 0.1" yr"' and agreeing 
within twice our uncertainties is 0.003. Therefore, we expect 
to find at most 0.0054 unrelated pairs similar to 2M1258AB 
in our entire search area, which implies a 99.5% probability 
of physical association for 2M1258AB. 

For thoroughness, we note that the angular separation of 
2M1258AB appears to increase monotonically over the three 
imaging epochs (2MASS, SDSS, CPAPIR) based on positions 
provided in table 1 . However, the trend falls within our uncer- 
tainties and is completely consistent with a constant separa- 
tion over time. Furthermore, upon examination of the 2MASS 
ATLAS images, we find that this apparent trend is broken 
when measuring the angular separation directly from the J 



and H band images, and only holds using the band image. 
Therefore, the trend appears to be coincidental, and it remains 
most likely that the system is bound. 

3. PHYSICAL PROPERTIES OF THE COMPONENTS 

A summary of observational and physical properties of 
2M1258AB is given in table|2] 

Spectral types for the components were determined primar- 
ily from the best-fitting referenc e spectra from the IRTF spec- 
tral library ICushing et al.ll20651) . Reasonable fits were found 
from M5 -M6.5 for the primary, and M6.5-M8 for the sec- 
ondary. Figure [T] shows our spectra of 2M1258AB plotted 
alongside those of other M4-M9 dwarfs for comparison. 

A weigh ted average o f the H 2O A, B, and C indices'* 
defined by iMcLean et al.l (12003 1) indicates spectral types of 
M6±0.8 and M7±0.8 for the primary and secondary respec- 
tively, in agreement with the former determination. Based on 
our fits to NIR reference spectra and the H2O spectral indices 
we assign spectral types of M6±l for the primary and M7±l 
for the secondary. The error bars are representative of the am- 
biguity present in assigning M dwarf spectral types based on 
low resolution NIR spectra, as is evident from figure[T] 

In order to determine absolute J, H, and, Ks magnitudes 
and their associated uncertaintie s we cornpiled a list of M 
dwarfs with measured pa rallaxes dJao et al.ll2005l; ICosta et al] 
120061: iHenrv et all 120061) . From this list, M dwarfs wifli 
spectral types within 0.5 subclass of the component spectral 
types were used to determine mean absolute magnitudes for 
2M1258AB, and the standard deviation of M dwarfs within 
1.0 subclass was used as a measure of the uncertainty. By 
comparing our derived absolute magnitudes to the 2MASS J, 

H, and Ks magnitudes we calculated spectroscopic distances 
for the components of 115 ±20 pc and 95 ± 18 pc for the 
primary and secondary respectively. This implies a projected 
separation of 6700 ± 800 AU at the average system distance 
of 105 ±13 pc. 

The equivalent widths (EWs) of the 1.25 fim K I and 

I. 20 /im FeH features can be used as rough gravity indica- 
tor s in NIR spec t ra. U sing the wavelength regions defined 
by iGorlova etall ( l2003l) we measured the K I EWs to be 
5.0 ±1.7, 10.6 ±1.6 and the FeH EWs to be 12.9 ±3.5, 
6.6 ±3.2 for the primary and t he secondary respect ively. 
When compared with figure 8 of IGorlova et al.l (1200 3h . the 
measured EWs-with the exception of the FeH EW in the 
secondary's spectrum-are most consistent with that of field 
dwarfs with \ogg > 4.5, and provide no evidence of youth. 

^ the M5 match was to Gl 866ABC, which has an alternate classification 
of M6 bv lReid et alj flOOj) 

■* we have chosen no to include the H2O D index as it predicts unreasonably 
early spectral types for both components, possibly a systematic effect of the 
telluric conection 
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TABLE 2 

Physical Properties of 2M 1 25 SAB 



Quantity 


A 


B 


2MASS Designation 


J12583501+4013083 


J12583798+4014017 


fia COS 5 (mas yr"' ) 


76±8 


76±8 


fJ-S (mas yr') 


-115±9 


-108 ± 10 


J (mag)" 


15.59 ±0.05 


15.61 ±0.05 


H (mag)" 


14.84 ±0.06 


14.89 ±0.06 


K, (mag)" 


14.43 ±0.06 


14.62 ±0.07 


r' (mag) ^ 


20.74 ±0.04 


22.00 ±0.07 


/' (mag) 


18.40 ±0.02 


19.18±0.03 


z' (mag) " 


17.18±0.02 


17.68 ±0.02 


Mj (magf 


10.20 ±0.70 


10.71 ±0.41 


Mh (magf 


9.55 ±0.64 


10.07 ±0.40 


MKs(mag)^ 


9.17±0.62 


9.70 ±0.37 


d(pc) 


115 ±20 


95± 18 


Spectral type 


M6±l 


M7±l 


T,ff (K)'J 


2850 ±300 


2620 ± 170 


Mass (Mq) 


105+" 


091+" '" 
"■"^'-0.007 


Angular separation (") 


63.38 ±0.05 


Projected separation (AU) 


6700 ±800 


Binding energy (10*' erg) 


0.25 ±0.08 


Mass ratio 


0.87 ±0.25 



" magnitudes from the 2MASS Point Source Catalog'' PSF magni- 
tudes from the SDSS DR6 Catalog'^ Estimates based on M dwarfs 
within ±1 spectral type with measured parallaxes'' based on the Teff 
versus spectral type relationship of Golimowski et al. (2004) 

However, the large errors due to the low S/N of our spectra 
make it impossible to draw firm conclusions about the age or 
gravity of 2M1258AB. Assuming an age of 1-5 Gyr (as an 
upper limit, see ® the comp onent masses can be estimated 
from evolutionary models (B araffe et alj|1998l) . Masses of 

+0.029 



0.105-;-^ Me 



and 0.09 lilio Mq were found for the pri- 



mary and secondary respectively by minimizing the error- 
weighted sum of square deviations between our derived abso- 
lute JHKs magnitudes and the models, with errors bars repre- 
senting 1(7 deviations from the minimum. This determination 
is also consistent with masses derived from effective tempera- 
tures rather than magnitudes. The mass ratio of 2M1258AB is 
0.87 ± 0.25, consistent with other wid e VLM binaries which 
tend to have mass ratios close to unity (iBurgasser et al.l2007h . 

4. DISCUSSION 

Of the previously known VLM binaries, only 2 have sepa- 
rations^ OOOAU: K61 AB with a separation of - 1800 AU 
(ICaballeroll2007bl) and 2M0126AB with a separation of - 
5100 AU dArtigau et arj|2007l) . In addition, a probable VLM 
binary with a separatio n of ~ 1700 AU has been identified by 
ICaballero et al.l (l2006h in the a Orionis cluster. With a sepa- 
ration of - 6700 AU 2M1258AB is the widest VLM binary 
yet discovered (see figure |2]). The existence of 3-4 of such 
ultrawide systems is strong evidence that they are not statis- 
tically rare chance alignments, but rather natural outcomes of 
star formation. 

Many other authors have discussed the implication s of wide 
VLM binaries for variou s formation models (e.g. iLuhmanI 
l2004t lArtigauet all 120071) . As with other wide VLM bina- 
ries, 2M1258AB would not have survived strong dynami- 
cal interactions during the f ormation process and is incon- 
sistent with ejection models jReipurth & Clark3l2001h . For- 
mation of the secondary in the disk of the primary (e.g. 
IStamatellos et al .'2007: Boss 2000) can also be strongly ruled 
out since disks around low ma ss stars tend to be at m ost a 
few hundred AU in diameter dVicente & AlvesI l2005h . and 




1.5 2.0 
Wavelength (^im) 

Fig. 1 .— Spectra for 2M1258AB taken with the SpeX prism at the IRTF, 
alongside compariso n spectra of field dwarfs from the IRTF spectral Ubrary 
ICushing et alj2005l) . 

tend to b e only a few perce nt as massive as the stars they 
surround (Schol z et al.l 120 06'). Thus this mechanism is un- 
likely to produce near-equal mass binaries, nor binaries with 
such wide separations. Alternatively, since it is not uncom- 
mon fo£jiighe£jiiass_st^_to form weakly bound binaries 
(e.g. lDuquennoy & Mayor! 1 99 Ih . models in which VLM stars 
form similarly to their more massive counterparts are better 
able to account for the ex istence of wide VLM pairs. Tur- 
bulent fragmentation (Pa doan & Nordlundll2002h is one such 
model, capable of producing low-mass cores down to ~ 3 
Jupiter masses. 

Higher order multiplicity may result in higher total system 
masses for wide VLM binaries, implying larger Jeans masses 
in formation. For more massive binaries, a large fraction 
25%) of wide (> 1000 AU) binaries are higher order multiples 
(iMakarov et al.ll2008h . Of the two previously known ultraw- 
ide VLM binaries, at least one, K61 AB, is a candidate triple 
system, where the secondary is a suspected spectroscopic 
double with total mass of 0.08 Mq dBasri & Reinersll2006t) . 
Other known examples among low-mass and VLM stars in- 
clude: LP 213-68A(BC) with M6.5 and M8 components an d 
a separation of 230 AU ( Gizis et al.ll200tf , 'Close et al."2003'): 
G 124-62ABab, M4.5/L1, 1900 AuTSeifahrt et al. 2005J; 
GJ 1 245ABC,(M5.5/M5.5)/M8, 38 AU (McCarthv_etd] 
IT988h: USco J160611.9-19 3532 AB, M5/M5,160b AU 
(Kra us & Hillenbrand! l2007bh : and c andidate triple LP 714 - 
37ABC M5.5/(M7/M8), 33 AU (Ph an-Bao et aH l2006l) . 
Based on 2MASS and SDSS photometry there is some evi- 
dence that 2M1258AB may also be a triple system. Although 
the primary clearly has an earlier spectral type than the sec- 
ondary (see figure[T]i, the apparent J, H, and K magnitudes of 
the components are nearly indentical. Referring to table 3 of 
iHawley et al.l (l2002h . we note that the difference in mean ab- 
solute J magnitude between M6 and M7 dwarfs is 0.74 mag, 
suggesting that the secondary may itself be an unresolved 
equal-mass binary. If both Kol AB and 2M1258AB turn out to 
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Fig. 2. — Sepai'ation and binding energy versus total system mass for known binary systems. Stellar binaries (dots) are from lTokovininH 1997|):|Fischer & Marcyl 
(1992); Reid & Gizis 1 1997); Close et al. ( 1990); Kiaus & Hillenbrand l2007b). VLM binaries (open squares, open circles) are from Burgasser et al. 120075; 
[Caballero et al. ( 2006); Luhman et al. ( 2009) and the VLM binary archive, maintained by Nick Siegler at http.V/www. vlmbinaries. org , Open circles identify VLM 
binaries that are known members of young associations. 



be triples then this would suggest that higher order multiplic- 
ity is common among the widest VLM systems. The second 
widest VLM binary, 2M0126AB, has yet to be checked for 
additional companions. 

Wide weakly bound binaries such as 2M1258AB are sub- 
ject to orbital evolution and eventual disruption as they travel 
through the Galaxy, due to encounters with stars and giant 
molecular clouds. Orbital evolution is an important consider- 
ation when assessing the probability that 2M1258AB formed 
in such a wide orbit, and also provides a rough constraint 
on the system age. The theoretical framework and numeri- 
cal results for the evolution and lifetimes of wide binary sys- 
tems in the solar neighborhood is provided by [Weinberg et alj 
([198^ . and can be easily scaled to apply to VLM systems 

i e.g.,|A rtigau et al. 2007). From figure la of Weinberg et al] 
19871) . the net effect of diffusive collisions is that an en- 
semble of wide systems will retain approximately the same 
average separation over time, provided they have not been 
disrupted. The dispersion in separation after ~ 1 Gyr is on 
the order of A log a « log 2. Thus the orbital separation of 
wide VLM binaries remains approximately within a factor of 
2 of the original separation. Even allowing for the possibility 
that binaries such as 2M1258AB as well as 2M0216AB were 
formed at half their current separation, they would still have 
initial separations an order of magnitude greater than that of 
other wide VLM binaries and present the same challenges for 
formation scenarios. The disruption lifetimes of wide bina- 



ries are also of interest as they can provide upper limits on the 
system age. An order of magnitud e estimate of the survival 
half-life, using figure 2 of Weinber g et al.l ( Il987l) for a binary 
with a/M,„, = 0.16 pc Mq~^ , yields fi/2 « 1 Gyr. If we as- 
sume the physical separation of the system to be greater than 
the projected separation then the survival half-life is slightly 
lower than this. Since it is unlikely that 2M1258AB is more 
than a few times older than fi/2, the upper limit on the system 
age of 1-5 Gyr used for deriving masses in ^is reasonable. 

Further characterization of this system would benefit from 
optical spectra to better constrain the spectral types and activ- 
ity of the components, as well as a parallax measurement to 
obtain a more accurate distance and hence separation. Both 
the primary and secondary should be checked for additional 
companions. 
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